Freezing and Collapse of Flexible Polymers 
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' We analyze the freezing and collapse transition of a simple model for flexible polymer chains on 

simple cubic and face-centered cubic lattices by means of sophisticated chain-growth methods. 
In contrast to bond-fluctuation polymer models in certain parameter ranges, where these two 
conformational transitions were found to merge in the thermodynamic limit, we conclude from 
our results that the two transitions remain well- separated in the limit of infinite chain lengths. 
The reason for this qualitatively distinct behavior is presumably due to the ultrashort attractive 
' interaction range in the lattice models considered here. 



1 Introduction 



a 

It is well known that single homopolymer chains undergo a structural coil-globule tonsi- 
ls) | tion at the so-called 0-point. Much theoretical, experimental, and algorithmic work was, 
and still is, spent to localize that point for various homopolymer models. In principle, 
there is no longer any difficulty to investigate lattice models up to very long chain lengths. 
■ A well working technique to deal with the problem is the (n)PERM algorithm. 1 ' 2 

Since a relatively short time, thanks to generalized-ensemble versions of PERM, 3 ' 4 it is 
also possible to investigate the temperature range far below the 6-temperature, where an- 
other transition called ground-state-globule, liquid-solid (crystallization) or freezing tran- 
sition can occur. 3 ' 5 A recent study of the bond-fluctuation model with respect to these 
different transitions for example showed, that the crystallization and the coil-globule tran- 
sitions may, but generally do not, coincide in the thermodynamic limit, depending on the 
p->) | interaction range. 5 

In this work, motivated by above mentioned studies, we will report results for the 
structural transitions of Interacting Self- Avoiding Walks (ISAW) on the simple cubic (sc) 
and face-centered cubic (fee) lattice, the simplest model for flexible, interacting polymers. 6 



2 Model and Methods 



In the ISAW model, a local attractive interaction between non-bonded nearest neighbors 
is assumed. The total energy of an ISAW is given by E — —n, where n is the number 
of such contacts. To simulate the model, we use the nPERM algorithm for the study of 
the 9-transition of very long chains and a generalized-ensemble version for the simulation 
of the low-temperature behavior of short polymers. The PERM method is a chain-growth 
algorithm based on the Rosenbluth method. 7 It includes population control by pruning or 
enriching populations during the growth, depending on threshold weights. 1 ^ 
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3 Results and Discussion 

For investigating structural transitions of the polymer model, we calculate the specific heat 
and analyze its peak structure. It is expected that even for polymers of finite length, peaks 
of fluctuating quantities signalize conformational activity. 

Figure QJleft) shows specific-heat peaks of polymers on the sc lattice with lengths 
8 < N < 125. We see that there is no uniform scaling behavior of the peaks as it was 
found for the bond-fluctuation model. 8 A first view does not show any regularity at all 
regarding the low-temperature peaks. However, the dependence of the peak temperatures 
on the chain length exhibits more systematics, see Fig.QJright). 
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Figure 1. Left: Map of specific-heat maxima for several chain lengths. Circles (0) symbolize the peaks (if 
any) identified as signals of the collapse (Tpmax > 1). The low-temperature peaks (+) belong to the excita- 
tion/freezing transitions (T^max < 0.8). Right: Collapse and excitation/freezing peak temperatures of the same 
specific-heat peaks. 

The freezing-transition temperatures show a sawtooth-like behavior which is due to opti- 
mal monomer alignment to the underlying lattice. At the lowest peak temperatures, we find 
chains with very compact ground states which are arranged as cubes or compact cuboids, 
respectively, e.g., for N — 27, 36, 48, etc. The corresponding chains have an energy gap of 
AE — 2 between the ground state and the first excited state. For the sc lattice, this can eas- 
ily be explained. The first excited state can be constructed by removing a monomer from 
the corner of the compact state (breaking 3 contacts) and placing it somewhere at the sur- 
face (gaining 1 contact). All these chains have a very pronounced freezing-transition peak 
(not shown). On the other hand, at the other side of the "teeth", the respective chains with 
one more monomer reside. Here, the formerly pronounced low-temperature peak becomes 
very weak. 

Let us briefly say some words on the 9-point. Despite of very precise measurements, 
the nature of this transition is not yet completely understood, considering for example 
predicted logarithmic corrections 9 which could not be resolved so far in numerical data. 
Figure [2] shows an illustrating picture of the studied objects (left) and data of transition 
temperatures for different chain lengths as well as fits to the data (right). There are several 
approaches to extrapolate T& by fitting T C (N). We used here the mean-field like fit 6 

1 ' 1 0) 



02 

N ' 



T C (N) T 

which was found to be consistent with numerical data obtained in grandcanonical analyses 
of lattice homopolymers and the bond-fluctuation model. 10 
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Figure 2. Left: Two typical conformations of a 4096mer on a sc lattice at different temperatures T -C 7©; 
E = —5367 and T 3> Tq; E = —1054. Right: Inverse collapse temperatures for several chain lengths on sc 
(N < 32 000) and fee lattices (N < 4 000). Drawn lines are fits according to Eq. {T}. Error bars are shown, but 
may not be visible. 



Optimal fit parameters using the data in the intervals 200 < N < 32 000 (sc) and 
100 < N < 4 000 (fee) were found to be T^ c = 3.72(1), a x » 2.5, and a 2 « 8.0 (sc), and 
Tq = 8.18(2), a± » 1.0, and a 2 ~ 5.5 (fee). This agrees very well with data published 
so far, 1 ' 11 but does unfortunately not solve the problem of correction terms either. A very 
detailed analysis including different fit ansatze and parameter estimates can be found in 
Ref. 6. 

To summarize, there exists a clear low-temperature freezing transition below the 
0-point, which is strongly influenced by lattice restrictions. It is shifted with increasing 
chain length to lower temperatures (Tq^(N) w 0.4) and jumps at chains with very com- 
pact ground states back to a value of T(yw(N) « 0.6. The temperature interval, in which 
the freezing peaks fluctuate, does not change when varying the chain length, whereas the 
finite-length 0-temperature does. 
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